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Abstract 

The primary present-day observables upon which theories of galaxy evolution are based are a system's 
morphology, dynamics, colour, and chemistry. Individually, each provides an important constraint to any 
given model; in concert, the four represent a fundamental (intractable) boundary condition for chemody- 
namical simulations. We review the current state-of-the-art semi-analytical and chemodynamical models for 
the Milky Way, emphasising the strengths and weaknesses of both approaches. 

Keywords: Galaxy: abundances — Galaxy: evolution — Galaxy: formation — Galaxy: kine- 
matics and dynamics — methods: n-body simulations 



1 Introduction 

Stars and interstellar gas in galaxies exhibit diverse chem- 
ical element abundance patterns that are shaped by their 
environment and formation histories. The aim of Galactic 
Chemical Evolution (GCE) is to use the observed abun- 
dances to unlock earlier epochs in the Universe, probe the 
mechanisms of galaxy formation, and gain insight into the 
evolution of stellar systems. 

Models for the chemical evolution of galaxies need to 
account for the collapse of gas and metals into stars, the 
synthesis of new elements within these stars, and the sub- 
sequent release of metal-enriched gas as stars lose mass and 
die. An additional feature of most models is the ongoing 
accretion of gas from outside the system. The most sophis- 
ticated models also incorporate a self-consistent treatment 
of the system's dynamics - both collisionless and dissipative 
components - either under idealised (semi-cosmological) con- 
ditions or within a full cosmological framework. Coupling 
GCE codes to a spectrophotometric evolution package fur- 
ther ensures that the models are constantly tested against 
real-world observational constraints. 

Semi-analytic homogeneous models make simplifying as- 
sumptions that enable the mean trends of galactic systems 
to be calculated by numerically solving a set of equations 
governing the formation, destruction, and distribution of 
the elements as they cycle through gas and stars. One 
strength of these models is that they typically have the 
fewest number of free-parameters, making convergence to 
a unique solution more likely. A weakness of homogeneous 
models is the inherent assumption that stellar ejecta from 
dying stars is instantly mixed back into the ambient inter- 
stellar medium (ISM). Inhomogeneous GCE models relax 
this so-called "instantaneous mixing approximation" in a 
semi-analytical manner, allowing consideration of observed 
trends in dispersion in various galactic observables. The 
self-consistent treatment of not only GCE, but the dy- 
namics of a galaxy's gas, stars, and dark matter, remains 
the purview of chemodynamical codes. Each of the above 
are complementary tools for deconstructing the formation 
and evolution of systems such as our Milky Way: semi- 
analytical models can cover a range of parameter space that 
a chemodynamical code cannot, due to the many orders of 
magnitude difference in the respective computational de- 
mands, while the latter afford a coupling of the dynamics 



of the system to that of the GCE, in a manner not otherwise 
available. 

In the review which follows, we present a biased overview 
of contemporary research in the field of Galactic Chemical 
Evolution. For seminal reviews tracing the development 
and principles of this topic, the reader is referred to Tinsley 
(1980), Matteucci (2001), and references therein. Section 2 
summarises the most popular formation and evolutionary 
scenarios and describes the relationships between different 
components of the Galaxy. Principles of homogeneous and 
inhomogeneous semi-analytical models are presented in § 3 
and § 4, respectively, while in § 5, the state-of-the-art in 
three-dimensional cosmological chemodynamical codes is 
reviewed. Coupling the GCE predictions from both of these 
approaches to galaxy evolution to the colour and luminosity 
information provided by spectrophotometric codes is dis- 
cussed in § 6. Potential future areas of interest to the field 
are itemised in § 7. 



2 Structure and Formation of the Galaxy 

The thin disk and bright inner bulge are the brightest com- 
ponents of the Milky Way. They are housed within a much 
more extended and diffuse spheroidal stellar halo - as dis- 
tinct from the dark matter halo. In addition, a fourth main 
component - the thick disk - was identified by Gilmore & 
Reid (1983). The latter found that star counts obtained 
by an extensive UK Schmidt photometric survey were well- 
fitted by a multi-component stellar population model con- 
sisting of 1) a thin disk, 2) thick disk, and 3) halo. This 
three-component model breaks down within the central kpc 
of the Galaxy where a dense metal-rich bulge begins to 
dominate (Gilmore, Wyse & Kuijken 1989). The thick disk 
component has an exponential scale height of 1350 pc - 
about four times greater than the scale height of the thin 
disk - and comprises ~2% of the nearby stars. Photomet- 
ric studies of external galaxies had already established that 
thick disks are common to spirals (van der Kruit & Searle 
1982). A popular explanation for the presence of thick disks 
is that mergers with smaller satellites during early times 
heated the thin disk (Wyse & Gilmore 1993). 

In order to distinguish individual populations of stars, 
one wishes to know not just their spatial distribution, but 
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Table 1: Milky Way Properties 
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Scalchcight Scalclcngth 




Age (Gyr) a 


[Fe/H]° 


(kpc) (kpc) 


Halo 


14 


-1.78 


Effective Radius~2.7 b 


Thick disk 


11 


-0.78 


~ 0.75 c 3.5 C 


Thin disk 


5-7 


-0.14 


~ 0.33 c 2.25 c 


Bulge 


10 





Effective Radius~1.2 d 



"Robin et al. (2003) 

fa de Vaucouleurs profile - Buscr et al. (1998) 
c Chen et al. (2001) 

d de Vaucouleurs profile - Yoshii & Rodgers (1989) 



their kinematics, chemical abundances, and ages. A com- 
plete dataset of this information should enable one to re- 
construct the formation and evolution of the Milky Way. 
A wealth of past, present, and future surveys and satellite 
missions (e.g. HIPPARCOS 1 , RAVE, 2 GAIA 3 ) offer (or 
will offer) the opportunity to determine the order in which 
the Galactic components formed, whether they evolved in- 
dependently of one another, and how important merging 
has been in assembling the Milky Way. The mean age and 
metallicity of the halo, thick disk, thin disk, and bulge are 
shown in Table 1. 

Traditionally two scenarios have competed to explain 
the formation of the Milky Way: 

1) The first scenario, proposed by Eggen, Lynden-Bell & 

Sandage (1962), describes the rapid monolithic col- 
lapse of a protogalactic gas cloud to form the halo. 
The Galactic disk would have subsequently formed as 
the residual gas disipationally collapsed. This would 
naturally give rise to two population of stars: an 
older, more metal-poor group found in the halo; and 
a younger, more metal-rich group orbiting closer to 
the Galactic mid-plane. 

2) Searle & Zinn (1978) offered an alternative to the mono- 

lithic collapse picture, proposing that the Galaxy was 
constructed from smaller cloud fragments, in which 
stars may have already started forming. 

The Galaxy's true formation history is likely to lie some- 
where between the two extremes of primordial collapse and 
hierarchical formation. Chemical properties of stars pro- 
vide important clues into disentangling the puzzle of the 
Galaxy's formation. The relative abundances of certain el- 
emental species act as "cosmic clocks", by which the for- 
mation timescales of various stellar populations can be de- 
termined. A popular cosmic clock is the ratio of an element 
like oxygen, that is born mostly in massive very short-lived 
stars, and an element like iron, whose creation is linked 
to lower mass longer-lived progenitors (Gilmore, Wyse & 
Kuijken 1989). 



3 Homogeneous Models 

Homogeneous GCE models have traditionally formed the 
cornerstone of this field and consequently have a rich liter- 
ature to draw upon. The basic ingredients, observational 
constraints, and several weaknesses are highlighted below, 

1 http : //astro . estec . esa . nl/Hipparcos/ 
2 http : //astronomy . swin . edu . au/RAVE/ 
3 http : //astro . estec . esa . nl/GAIA/ 



although it should be emphasised that much of this dis- 
cussion pertains also to the inhomogeneous and chemody- 
namical models described in § 4 and § 5 (as many of the 
ingredients are common to all techniques) . 

3.1 Basic Ingredients 

The main ingredients of homogeneous GCE models are out- 
lined below and we discuss their relationship to one another 
through the basic set of chemical evolution equations. 

• Stellar Yields and Lifetimes: Almost all elements 
heavier than helium originate from stars. Stars enrich 
the ISM with their own unique pattern of elements 
depending on their mass and initial metallicity. The 
predicted stellar yields consequently form the back- 
bone of the study of GCE. The grids of yields utilised 
in the models are the outcome of computations of 
stellar evolution and vast networks of nuclear reac- 
tions. For the purposes of chemical evolution, stars 
are often divided into three categories: 

Massive stars (m > 10 M Q ) evolve quickly because 
their enormous gravitational potential acceler- 
ates the nucleosynthesis process. Their death 
is marked by a violent supernova (SN), leaving 
behind a neutron star or black hole. Although 
massive stars are much rarer than their lower 
mass counterparts, they are the main source 
of most of the heavy elements (i.e. metals) 
in the Galaxy. Figure 1 shows the production 
factors from massive stars predicted by the de- 
tailed nucleosynthesis calculations of Woosley 
& Weaver (1995, left panel) and the FRANEC 
code (kindly provided by A. Chieffi 2003, priv 
comm; FRANEC is described in Chieffi & Limongi 
2002, right panel). 4 The dotted line at [X/O] =0 
indicates the solar elemental abundance pattern 
relative to O. The solar abundance pattern of 
most metals is adequately reproduced by mas- 
sive stars but C, N and the iron-peak elements 
require additional production sites. 

Intermediate- and low-mass stars (ILMS) live longer 
than their massive star counterparts, due to 
their lower density. They greatly outnumber 
the heavier stars but do not produce signifi- 
cant quantities of many elements besides he- 
lium, carbon, nitrogen, and certain isotopes, 
which are created through hydrostatic burning 
and expelled in stellar winds and planetary neb- 
ulae (e.g. van den Hoek & Groenewegen 1997). 
Very low mass stars (m < 1 M©) have lifetimes 
comparable to the age of the Galaxy and there- 
fore serve to lock up the gas supply. 

ILMS in binary systems may culminate in powerful 
supernovae explosions classed as Type la. The 
exact physical mechanisms behind SNe la are 
still an open question but one popular theory 
holds that the mass lost by a binary star as it 
evolves is accreted by its smaller white dwarf 

4 A generic comparison of massive star yields can be found in 
Gibson, Locwenstein & Mushotzky (1997). 
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(WD) companion until the WD can no longer 
be sustained by electron degenerate pressure. 
Then the entire mass of the WD is ejected in 
a violent explosion that converts much of the 
stellar material into iron (e.g. Iwamoto et al. 
1999). 

Initial Mass Function (IMF): The precise form 
of the IMF dictates the number of stars born in a 
given mass interval in each generation of stars. This 
in turn sets the rate at which different elements are 
released into the ISM, thus influencing both the rela- 
tive and absolute elemental abundances. Most IMFs 
in the literature consist of simple single- or multi- 
component power law specified over a mass range 
from m ~ 0.1 M to an upper mass limit typically 
between 40-100 M Q (e.g. Salpeter 1955; Scalo 1986; 
Kroupa, Tout & Gilmore 1993). 

Star Formation Rate (SFR): While laws of star 
formation can be calculated from first principles, chem- 
ical evolution models invariably use a functional form 
that has been derived empirically. An ample supply 
of gas is the first condition needed for star forma- 
tion, so it is not surprising that one of the simplest 
laws has SFR oc <jg as , where o- gas is the surface den- 
sity of gas and the exponent k may range from 1-2 
(Schmidt 1959). Other star formation laws presume 
that factors such as total mass density and/or Galac- 
tocentric radius play a role. Dopita & Ryder (1994) 
found that a law given by SFR oc o-% 3 s a^f al satisfac- 
torily describes the correlation between Ha emission 
and I-band surface brightness in spiral galaxies. 

Gas Flows: In the simplest scenario, our model 
Galaxy in each radial bin can be considered a closed- 
box consisting of primordial gas from which stars are 
born according to the chosen SF and IMF prescrip- 
tions. In the classic closed-box model (e.g. Pagel 
& Patchett 1975) there is no gas loss or gain; at 
time = all Galactic matter is present as primordial 
gas from which stars immediately form. This type 
of model, characterised by an intense period of early 
star formation, provides a reasonable account of the 
formation of the halo and bulge of the Milky Way. 
When applied to the Galactic disk, however, the ba- 
sic closed-box model leads to an excess of metal-poor 
stars with respect to the observed metallicity distri- 
bution of nearby long-lived stars: the so-called "G- 
dwarf problem" (Pagel & Patchett 1975). More re- 
alistic models overcome this problem by allowing the 
Galactic disk to form via continual accretion of gas. 
The infall rate as a fucntion of radius r and time t 
often takes the form: 

j t o gas {r,t) = A{r)e- t '^\ (1) 

where r(r) is the exponential infall timescale and 
the coefficient A(r) must satisfy the constraint that 
f now A(r)e~ t/T( - r ^ = o-gas(r), where a 9as (r) is the 
present-day surface density profile. Many elliptical 
and dwarf galaxies are expected to have undergone 
periods in which the energy from supernova explo- 
sions exceeds the gravitational energy of the system, 



causing the interstellar gas to be expelled in a galac- 
tic wind (Larson 1974; Dekel & Silk 1986; Matteucci 
& Tornambe 1987; Gibson 1997; Recchi, Matteucci, 
& D'Ercole 2001). However, the outflow of gas is 
not expected to feature in the history of most spiral 
galaxies and is usually neglected in the models. Ra- 
dial flows of gas have also been investigated within 
the homogeneous GCE paradigm (e.g. Portinari & 
Chiosi 2000). 

• Galactic Components: The infall of gas described 
above might occur in multiple episodes that corre- 
spond to the formation of the individual halo, bulge, 
thick disk and thin disk components. Recent works 
have adopted different premises regarding the extent 
to which these components interact with each other 
through gas exchange. For instance, enriched halo 
gas might be funneled into the bulge or else might 
collapse further to form the disk. Alternatively, the 
halo and thin disk might evolve coevally yet inde- 
pendently. The answer to these questions is partly 
clouded by challenges in distinguishing the various 
stellar populations. The properties that we wish to 
measure in each population, such as metallicity and 
kinematics, are also the properties used to define the 
populations. 



The chemical composition of the interstellar medium as a 
function of time and Galactocentric radius is described by 
equations that balance: 1) processes that deplete chem- 
ical elements from the interstellar medium, namely star 
formation and perhaps Galactic winds; and 2) processes 
that replenish the ISM, such as stellar winds, SNe, and in- 
falling gas. The GCE equations can be solved analytically 
given the assumption that stars release their ejecta instan- 
taneously at the time of their birth. This approximation 
is reasonable only for very massive and short-lived stars 
and it precludes one from reproducing the relative trends 
of elements that are restored to the ISM by different mass 
stars on different characteristic timescales. For GCE mod- 
els to have the power to interpret the variation of diagnostic 
abundance ratios, the "instantaneous recycling approxima- 
tion" must be relaxed. There is another assumption, how- 
ever, that few GCE models relax: namely, that the stellar 
yields are instantaneously mixed into the ISM. For this rea- 
son, homogeneous GCE models are not the ideal tools for 
studying the Galaxy in its infancy, when chemical abun- 
dances in the ISM could be traced to individual supernovae 
(or for present-day observables that originate from that era 
- e.g. the Galactic halo). We return to this assumption in 
§ 4. 

Standard models of spiral galaxy chemical evolution as- 
sume azimuthal symmetry and collapse the object into a flat 
disk such that radius is the only spatial variable. Defining 
<Ji(r,t) as the mass surface density of species i at radius r 
and time t, then the rate of change of Oi (r, t) is given by: 



dt 



<Ti{r,t) = I ip(r,t-T m )Yi(m,Zt- T „ 



)^M dm 



+ ^0"i( r > t)infall 



3.1.1 The Equations 
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Figure 1: Production factors relative to O on a solar logarithmic scale from a single generation of massive 
stars using the metallicity-dependent yields of Woosley & Weaver 1995 (left panel) and those of FRANEC 2003 
(right panel). The latter were kindly provided by A. Chieffi (2003, priv comm). Yields were integrated over a 
Salpeter (1955) IMF from 12-40 M©. The dashed line indicates the solar values (where log(N /N H )o + 12 = 
8.73 (Holweger 2001)) and dotted lines indicate deviations from scaled solar by a factor of two. For both sets 
of yields C, N, and some of the iron-peak elements are subsolar because they require additional sources such as 
lower mass stars and Type la SNe. The strength of the "odd-even" effect increases with decreasing metallicity 
in both cases, however the effect is more pronounced for FRANEC 2003. 



-X t (r,t)i>(r,t) 

— -^-Oi(r,t) ou t flow (2) 

where the four terms on the right-hand side of equation (2) 
correspond to the stellar ejecta, gas infall, star formation, 
and gas outflow rates, respectively, ip is the SFR, Yi(m, Zt- Tm ) 
is the stellar yield of i (in mass units) from a star of mass 
m and metallicity Z t - Tm , 4>( m ) is the IMF by mass, and 
Xi is the mass fraction of element i. By definition, the sum 
of Xi over all i is unity, and the total surface mass density 
is identical to the integral over the infall and outflow rates. 
rrii ow and m up are the lower and upper stellar mass limits, 
respectively, and r m is the lifetime of a star of mass m. 
In practice, the first term is split into three equations that 
deal separately with ILMS, Type la SNe progenitors, and 
massive stars (see Greggio & Renzini 1983 and Matteucci 
& Greggio 1986 for details). 

3.2 Observational Constraints 

The most thoroughly observed and best understood galaxy 
is the Milky Way, and in particular, the "local" solar neigh- 
bourhood. So extensive is the Milky Way dataset, many 
model ingredients can be well-constrained empirically. Thus 
our own Galaxy is often the gauge by which chemical evo- 
lution models are calibrated. Indeed, studies of the cosmic 
evolution of disk galaxies often adopt scaling laws based on 
the Milky Way (e.g. Boissier & Prantzos 2000). The ex- 
tent to which the Milky Way's IMF, star formation law, and 
nucleosynthetic behaviour can be applied to other types of 
objects at earlier epochs depends on how universal these 
prescriptions are and on whether we live in a prototypical 
galaxy. A minimal set of observational constraints for GCE 
models is described below. 



• Solar abundance pattern: Any chemical evolution 
model should be able to reproduce the solar abun- 
dance pattern, i.e. the pattern in the ISM 4.5 Gyr ago 
at the radius where the Sun was born. The Sun is the 
single star with the most complete set of abundance 
measurements. For consistency, one can also compare 
the chemical composition of solar system meteorites 
with estimates based on stellar spectral lines. As 
discussed above, the predicted solar enrichment pat- 
tern is chiefly controlled by the yields released in 1) 
Type II SN explosions of massive stars, 2) planetary 
nebulae and stellar winds of ILMS, and 3) Type la 
SN explosions of binary systems of ILMS. In addition, 
the predictions are sensitive to the stellar mass distri- 
bution (i.e. the IMF) and the SFR, since these set the 
relative contribution and enrichment timescale from 
different types of stars. A seminal attempt to simu- 
late the evolution of all elements up to and including 
zinc was made by Timmes, Woosley & Weaver (1995, 
hereafter TWW95) using the comprehensive grid of 
mass- and metallicity-dependent yields calculated by 
Woosley & Weaver (1995). Impressively, TWW95 re- 
produced the abundances of most isotopes in the Sun 
to within a factor of two. 

It is worth noting that even our understanding of 
the solar abundance pattern still has the power to 
surprise us. As recently as a decade ago, the solar 
oxygen abundance was assumed to be log(No/NH)© 
+ 12 = 8.93 (Anders & Grevesse 1989); it has be- 
come clear recently though that accounting for so- 
lar granulation and non-LTE effects properly leads 
to a significantly radical downward revision in the 
Sun's oxygen abundance by almost a factor of two (to 
log(N o /N H ) + 12 = 8.69 - Allende Prieto, Lambert 
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& Asplund 2001). Such a shift partially resolves the 
long-standing dichotomy between the Sun's oxygen 
abundance and that of the local ISM. 

• G-dwarf distribution: A much stronger constraint 
on GCE models is the distribution of stars as a func- 
tion of metallicity, since this represents the convolu- 
tion of the age-metallicity relationship and the star 
formation history. In order to probe the early Galaxy, 
one needs a sample of low mass stars such as G- or 
K-dwarfs whose main-sequence lifetimes are compa- 
rable to, or older than, the age of the Universe. GCE 
models have demonstrated that the paucity of low- 
metallicity dwarf stars can be explained if the Galac- 
tic halo formed first on a rapid timescale, followed 
by a slow build-up of the thin disk (e.g. Chiappini 
et al. 1997; Alibes et al. 2001; Fenner & Gibson 
2003). The excess of metal-poor stars predicted by 
simple closed-box models can also be avoided by as- 
suming prompt initial enrichment, perhaps by a first 
generation of extremely massive Population III stars. 
If this were the case, then the abundance pattern of 
Pop III ejecta should be evident in the lowest metal- 
licity stars. It has also been suggested that there are 
no very metal-poor stars left because they all had rel- 
atively short lifetimes due to low metallicity environ- 
ments favouring the formation of higher mass stars 
(e.g. Nakamura & Umemura 2001). A further con- 
sideration is that an initially pristine zero-metallicity 
star might have had its surface layers polluted by 
the accretion of metals from the ISM over the past 
~12 Gyr (e.g. Shigeyama, Tsujimoto & Yoshii 2003). 

• Evolution of abundance ratios: If elements Xi 
and X2 have different origins and different character- 
istic timescales for release into the ISM, then [X1/X2] 
vs [X2/H] acts as a clock by which chemical evolu- 
tion can be measured (e.g. Wyse & Gilmore 1988). 
For example, readily observable features of oxygen 
and iron in stellar spectra has encouraged the wide 
use of [O/Fe] vs [Fe/H] to diagnose the overall star 
formation history of galactic systems. As with most 
heavy elements, oxygen is produced chiefly in mas- 
sive and short-lived stars. Thus, oxygen enrichment 
immediately follows the onset of star formation. In 
contrast, at least half of the iron in the Galaxy prob- 
ably originated from Type la SNe (e.g. Alibes et al. 
2001), whose lower-mass and longer-lived progenitors 
introduce a time delay for iron enrichment. The re- 
maining iron comes largely from Type II SNe. The 
combination of high [Fe/H] and high [O/Fe] is under- 
stood to arise in systems that formed stars so rapidly 
that high metallicities were reached before SNe la 
had a chance to lower the [O/Fe] value (Smecker- 
Hane & Wyse 1992). Similarly, one might interpret 
low [Fe/H] and low [O/Fe] as a sign of a slowly evolv- 
ing system. 

• SFR and SN rates: The present-day star forma- 
tion and Type II and Type la SNe rates must be 
matched by a successful chemical evolution model. 
However these are fairly weak constraints given that 
we can only be reasonably certain about the current 
SFR and mean past rate. Finer details of the Galac- 
tic star formation history are difficult to recover and 
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are quite uncertain. The most direct way to infer 
the SFH is by determining the age distribution of 
stars; a method that relies upon unreliable stellar 
ages and assumptions about the IMF, stellar evolu- 
tion, scale height corrections, and stellar kinematics 
(Rocha-Pinto et al. 2000a). This technique is also 
somewhat circular, in that a SFH must have been 
assumed in order to derive the IMF. 

• Age-metallicity relationship: This is an impor- 
tant constraint, but again, a weak one given that 
the scatter in the observations (e.g. Ibukiyama & 
Arimoto 2002) can accommodate most model pre- 
dictions. Moreover, the very existence of an AMR, 
which had been well established by earlier studies (e.g 
Twarog 1980; Edvardsson et al. 1993; Rocha-Pinto 
et al. 2000b), has recently been challenged by inves- 
tigations demonstrating large intrinsic scatter and no 
significant trend of metallicity with age (e.g. Feltz- 
ing, Holmberg & Hurley 2001). 

• Gas and Abundance gradients: It has long been 
known that the Milky Way is more metal-rich toward 
its centre and more metal-poor at large Galactocen- 
tric distance (Tinsley 1980, and references therein). 
Using the oxygen abundance observed in H II regions 
and OB stars to trace metallicity, a metallicity gra- 
dient of —0.07 dex kpc -1 has been established (e.g. 
Smartt & Rolleston 1997). The abundance of met- 
als in a region of gas is particularly sensitive to the 
balance between the star formation and gas accre- 
tion rates. Therefore the predicted metallicity gra- 
dient of the Galactic disk depends strongly on how 
the star formation prescription and gas infall rate 
are assumed to vary with radius. Good fits to the 
data are obtained by "inside-out" formation scenar- 
ios, whereby the innermost disk is built-up on the 
shortest timescale (e.g. Larson 1976; Chiappini, Mat- 
teucci & Gratton 1997). Portinari & Chiosi (1999) 
showed that a SFR such as the Schmidt (1959) law, 
which varies only with the gas surface density, pro- 
duces a radial abundance profile that is too flat unless 
one assumes an unreasonably large variation in for- 
mation timescale from the inner to outer disk such 
that the far disk would currently be accreting at 
much higher rates than observed. The Dopita & Ry- 
der (1994) law, with a mild dependence on total mass 
surface density, yields a better fit to the metallicity 
gradient (Portinari & Chiosi 1999). The theoretical 
metallicity distribution of long-lived stars as a func- 
tion of Galactocentric radius is shown in Figure 2, 
assuming a Dopita & Ryder (1994) star formation 
law. 

• Isotopic abundances: Traditionally, chemical evo- 
lution studies have been concerned with monitoring 
the total abundance (or dominant isotope) of spe- 
cific elements in order to unravel the Galaxy's history. 
Recent advances in instrumentation have paved the 
way for research into individual isotopes that pro- 
vides new challenges for nucleosynthesis theory. For 
instance, Type II SNe models appeared capable of ex- 
plaining the magnesium isotopic ratios in intermedi- 
ate to solar metallicity stars, but the results from the 
solar neighbourhood model shown in Figure 3 reveal 
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Figure 2: The metallicity distribution function (MDF) 
of long-lived stars predicted by a single infall model 
at three different Galactocentric radii; 4, 8, and 14 
kpc. Open squares with error bars correspond to the 
MDF of nearby K-dwarfs (i.e. at a radius ~ 8-8.5 kpc) 
observed by Kotoneva et al. (2002). The theoretical 
MDF has been convolved with a Gaussian of disper- 
sion a = 0.1 dex in [Fe/H] to simulate the empirical 
uncertainties. 




Figure 3: The predicted metallicity evolution of mag- 
nesium isotopic ratio 26 Mg/ 24 Mg for the solar neigh- 
bourhood both with (solid curve) and without (dot- 
ted curve) a contribution from intermediate-mass AGB 
stars (from Fenner et al. 2003). Stellar data are from 
Gay & Lambert (2000, blue circles) and from halo and 
thick disk stars from Yong et al. (2003, orange dia- 
monds). The red square denotes the solar ratio. 



an underproduction of 26 Mg/ 24 Mg at low metallicites 
(dotted line) with respect to the latest data. The 
missing piece of the puzzle may be low metallicity in- 
termediate mass stars on the asymptotic giant branch 
(AGB), whose helium-shells may be hot enough to 
generate 25 Mg and 26 Mg by triggering a-capture onto 
22 Ne (Karakas & Lattanzio 2003). As Fenner et al. 
(2003) have demonstrated for the first time, the data 
at [Fe/H] < — 1 is much better matched after incorpo- 
rating the Karakas & Lattanzio AGB nucleosynthe- 
sis calculations in a chemical evolution model (solid 
line). Such detections of isotopic ratios in field stars 
and globular clusters may reveal a great deal about 
the relative role of different types of stars in various 
environments. 

3.3 Uncertainties and Weaknesses 

• Iron-peak Yields: Iron-peak elements are buried 
deep within the cores of massive stars near the ra- 
dius that separates the ejected material from the rem- 
nant. The location of this so-called "mass cut" is a 
free-parameter in stellar models, one which controls 
the relative abundances of the iron-peak elements as 
well as the Xi/Fe ratio in the ejecta. Abundances in 
very metal-poor stars can help constrain the choice 
of mass cut, however in order to simultaneously eject 
iron-peak elements in the correct proportions and re- 
cover the high observed [a/Fe] ratios, models need to 
incorporate mixing and fallback (Umeda & Nomoto 
2002) or asymmetrical explosions. Multi-dimensional 
simulations of explosive nucleosynthesis may reveal 
more about these processes (Travaglio, Kifonidis & 



Miiller 2003). 

• Shape and evolution of the IMF: A time-invariant 
IMF remains the best choice for modelling the general 
evolution of our own Galaxy (Chiappini, Matteucci 
& Padoan 2000), but peculiar abundance patterns in 
extremely metal-poor stars support the notion that 
the first generation of stars was biased towards higher 
masses (Chieffi & Limongi 2002). The upper limit of 
the stellar IMF, m up , is also uncertain and impacts 
upon the total amount of metals produced by each 
stellar generation. Figure 4 illustrates the the sen- 
sitivity of metal-growth in the solar neighbourhood 
to the upper IMF mass. Increasing m up from 40 Mg 
to 100 Me is expected to raise the metallicity by as 
much as 30%. This is due to the steeply increasing 
yield of O (which is the most abundant metal in the 
interstellar gas) as a function of stellar mass. 

• Black hole mass limit: Essentially a free param- 
eter in Type II SNe models, the mass above which 
most of the stellar material collapses to form a black 
hole must be addressed a posteriori, in an empirical 
manner (to recover the observable constraints alluded 
to earlier). Black hole progenitors are expected to re- 
lease most of their oxygen, carbon, etc., via pre-SN 
stellar winds, but heavier elements such as iron are 
expected to fall back onto the remnant. One may 
expect the mass range corresponding to black hole 
collapse is sensitive to metallicity (e.g. Maeder 1992) . 

• Limited dataset: Many conclusions about cosmic 
chemical evolution have been drawn in part from 
studies of the Milky Way, but the Milky Way is only 
one object. How unique or typical the Milky Way is 
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Figure 4: Evolution of Type II SN rate (red lines) 
and mctallicity (blue lines) predicted by a single in- 
fall phase model of the solar neighbourhood for an 
IMF upper mass limit of 40 M Q (dotted lines) and 
100 M(T)(solid lines). Mctallicity is defined as the mass 
fraction of elements heavier than He in the interstellar 
medium (ISM). The SN II rate is barely affected by 
changes in the upper limit to the (Kroupa et al. 1993) 
IMF, however metallicity is about 30% higher for the 
m up =100 Mq case. This reflects the steep increase in 
the yield of O (which is the dominant element con- 
tributing to metallicity) as a function of initial stellar 
mass. 



amongst other galaxies is unknown and it can be dan- 
gerous to take the Milky Way's evolutionary path as 
representative of most spirals. For instance, a system 
such as M31, that resembles the Milky Way in terms 
of size and morphology, shows evidence for dramat- 
ically different metallicity distributions of its stellar 
populations (Worthey & Espana 2003). As an ex- 
ample of the perils of drawing conclusions based on 
the solar neighbourhood, the almost constant ratio 
of Zn/Fe versus metallicity might be taken as a sign 
that zinc, like iron, owes much of its production to 
Type la SNe (Matteucci et al. 1993; Mishenina et al. 
2002). It is possible, however, that zinc comes pri- 
marily from Type II SNe, with a strongly metallicity- 
dependent yield that mimics the time-delay associ- 
ated with a SN la source (e.g. Timmes et al. 1995). 
As a consequence, zinc would not always vary in lock- 
step with iron but would depend on specific forma- 
tion histories. It is worthwhile to bear this in mind, 
as zinc is often used as a proxy for iron in the high- 
redshift universe. 



4 Inhomogeneous Models 

The halo has special significance for the formation of the 
Galaxy. If halo stars form a collisionless system, their or- 
bits contain information about the dynamics at the time 



of the formation. In addition, the abundances observed 
in low mass stars reflect the chemical abundances and in- 
homogeneities during halo formation. Element abundance 
ratios as a function of metallicity (generally) show increas- 
ing scatter with decreasing metallicity, 5 whereas at higher 
metallicities the scatter decreases to reach a mean element 
abundance which corresponds to the ratio of the stellar 
yields integrated over the IMF (the regime in which the 
homogeneous GCE models of § 3 are most applicable). 

The enrichment of the halo mainly depends upon the 
number of SN explosions and the manner in which the 
ejected gas is mixed with the ambient ISM. If the mix- 
ing volume is sufficiently large, the enrichment could be 
spatially homogeneous. Conversely, should mixing be in- 
efficient, significant (localised) abundance inhomogeneities 
could exist; in theory, gas in the vicinity of a SN might even 
bear the chemical imprint of that single event. If the lat- 
ter were the case, second-generation, extremely metal-poor 
halo stars may show an abundance pattern which matches 
the nucleosynthetic yields of a single Pop III SN. An highly 
inhomogeneous halo would naturally result in this case. 
Simulating the temporal and spatial history of chemical in- 
homogeneities in the halo of the Milky Way, through semi- 
analytical GCE models, is an extremely active field at the 
present. Important earlier work in this field include investi- 
gations by Malinie, Hartmann & Mathews (1991), Malinie 
et al. (1993), Pilyugin & Edmunds (1996), Copi (1997), van 
den Hoek & de Jong (1997), Ikuta & Arimoto (1999), and 
Travaglio, Galli, & Burkert (2001). In what follows, we 
review three of the more recent approaches to modelling 
inhomogeneous chemical evolution. 

4.1 Argast et al. (2000,2002) 

Argast et al. (2000,2002) evolve an ~15 kpc 3 region of mass 
10 8 Mq with a spatial resolution of 50 pc. At each time step 
(1 Myr), a certain number of cells are randomly chosen, and 
each one can create a star with a probability proportional to 
&gas (akin to the Schmidt Law described in § 3.1), under the 
assumption of a Salpeter (1955) IMF. Type la SNe were not 
included in their models, and the simulations were halted 
once [Fe/H] reached 1/10 solar. Stars which form out of 
material enriched by a single SN will inherit its abundances 
and show an elemental pattern which reflects the progenitor 
mass. The Argast et al. (2000) model has minimal ISM 
mixing, as the expansion of the SN remnant is the only 
dynamical process taken into account. 

The results of their Inhomogeneous Galactic Chemical 
Evolution (IGCE) modeling - in relation to the trend of 
abundance ratio scatter as a function of [Fe/H] - are in good 
agreement with observations, except for [Cr/Fe], [Mn/Fe], 
and [Ni/Fe]. 6 Argast et al. conclude that for [Fe/H]< 
—3, the ISM is unmixed and dominated by local inhomo- 
geneities polluted by individual SN events. For [Fe/H]> 
— 2, their model halo ISM reflects a true IMF-averaged 
abundance pattern and is considered "well-mixed" . Argast 
et al. do note though that an individual SN event can still 

5 As Spite et al. (2003) note though, this canonical "wisdom" 
needs some revision. 

6 The scatter in chromium and manganese does not increase 
(observationally) with decreasing iron abundance, while the ob- 
served scatter in nickel exceeds that predicted by their models. 
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have an impact upon a well-mixed ISM, leading to finite 
dispersions in abundance patterns at disk-like metallicities. 

4.2 Oey (2000,2003) 

The initial condition of Oey's (2000) IGCE model is that of 
a metal-free closed-box, in which the first generation of star 
forming regions is randomly distributed, occupying a vol- 
ume filling factor Q. Each region behaves as a superbubble 
powered by its enclosed Type II SNe. Oey then considers n 
subsequent generations of star formation, allowing the star 
forming regions to overlap. The main conclusion drawn 
by Oey (2000) is that the evolutionary state of a system 
is characterized by the product nQ, with the relative fill- 
ing factor of contamination having the same importance as 
the number of contamining generations. The high metal- 
licity tail of the MDF may provide a useful discriminant 
between the classical "Simple Model" and the Simple Inho- 
mogeneous Model proposed by Oey. The latter agrees with 
both the Galactic halo and bulge MDFs by varying only 
this single parameter nQ, with Q and n independent and 
roughly associated with the global star formation efficiency 
and age. 

The Simple Inhomogeneous Model assumes no large 
scale mixing beyond the superbubble radii, with metals 
uniformly distributed within the volumes of the hot su- 
perbubbles and cooling locally. Once mixing is allowed, 
the metallicity would be reduced by dilution, and an in- 
crease in nQ would be required to attain a given present 
day metallicity. Oey (2003) takes into account interstel- 
lar mixing processes in ordinary multiphase ISM (mainly 
diffusion and turbulent mixing 7 ), which does lead to lower 
metallicities within the superbubble. As a result of this 
analysis, turbulent mixing would appear to be more effi- 
cient than diffusion, but the lowered metallicity for parent 
enrichment events requires more evolution (higher nQ) to 
match the observed metallicities, and this in turn implies 
that the system's MDF drops off too steeply to match the 
data in the high metallicity tail of the Galactic halo MDF. 
Furthermore, the Pop III stellar fraction is too high com- 
pared to the observations, suggesting that a discrepancy 
remains between the model and the observations. Future 
developments in the model are eagerly anticipated. 

4.3 Tsujimoto et al. (1999) 

Under the Tsujimoto et al. (1999) formalism, halo star for- 
mation is confined to separate clouds of mass M c . Each 
cloud is initially composed of Pop III stars, with mass frac- 
tion xin, and gas that has yet to form stars. Subsequent 
generations of stars are assumed to form in SN remnant 
(SNR) shells. The mass fraction of each shell that turns 
into stars (e) is assumed to be constant. Heavy elements 
ejected from a SN event are assumed to be trapped and 
well-mixed within the SNR shell. Some of these elements 
go into stars forming of the next generation. This process 
continues until remnants are no longer capable of sweep- 
ing up sufficient gas to form shells. The mass of a shell 
M s h(m,t) — M e j(m) + M BW (m,t), where M e j(m) is the 
mass of the SN ejecta, and M sw (m, t) is the mass swept 

7 Note though that mixing between discrete enrichment re- 
gions is not allowed in Oey's (2003) prescription. 



up by a shell (assumed to be 6.5xl0 4 Mq throughout their 
analysis) . 

The free parameters of the model are the mass fraction 
xin of metal free Pop III stars initially formed in each cloud 
and the mass fraction e of stars formed in each SNR. These 
values are chosen to reproduce the observed [Fe/H] distri- 
bution function of halo field stars for [Fe/H]< — 1. If xm is 
too large, the total gas swept up by the first SNRs exceeds 
the entire amount of available gas, and the star formation 
stops at the first or second generation. If the process is to 
continue as a sequence of SN-induced star formation, xm 
must be ^10 -2 and e confined to a narrow range such that 
of order one massive star is born from each SNR. If e is too 
high, star formation soon stops with little enrichment. If e 
is too low, star formation will proceed until most of the gas 
is used up, with an excess of enrichment. 

Tsujimoto et al. (1999) conclude that the probabil- 
ity pin of observing a Pop III star amongst the general 
background field of halo stars, under the assumption that 
M c = 10 6 - 1O 7 M , should be 10~ 3 - 10~ 4 - consistent 
with current observational limits (Beers 2000). Their IGCE 
model also naturally recovers the frequency distribution of 
stars in the [Eu/Fe] — [Fe/H] plane; [Eu/Fe] spans >2dex 
for [Fe/H]< — 2, converging to a plateau by [Fe/H]~— 1. 
Further enhancements and applications of the model are de- 
scribed in Suzuki & Yoshii (2001), Tsujimoto & Shigeyama 
(2002), and Tsujimoto et al. (2002). 



5 Chemodynamical Models 

GCE is intimately related to the Galactic star formation 
history, and star formation is equally linked to the dynam- 
ical evolution of the Galaxy. The self-consistent treatment 
of the chemical and dynamical evolution of a system has 
long been recognised as desirable, but the computation- 
ally intensive nature of the simulations made this desire 
difficult to realise. Advances in hardware and numerical 
methods over the past decade though has finally allowed 
chemodynamical codes to realise their theoretical promise 
(e.g. Larson 1976; Samland, Hensler & Theis 1997; Car- 
raro, Lia & Chiosi 1998; Nakasato & Nomoto 2003; Brook 
et al. 2003a). Recent three-dimensional chemodynamical 
evolution codes are now being applied routinely to study 
disk galaxies such as the Milky Way. For example, Stcin- 
metz & Muller (1995) succeeded in distinguishing the chem- 
ical properties between halo, bulge, and disk stars (see also 
Bekki & Chiba 2001). Raiteri, Villata & Navarro (1996) 
and Berczik (1999) took into account metal enrichment by 
both Type II and Type la SNe, reproducing the correla- 
tion between [O/Fe] and [Fe/H] for stars in the solar neigh- 
bourhood. Brook, Kawata & Gibson (2003b) studied the 
metal enrichment from intermediate mass stars and dis- 
cussed how the distribution of carbon, nitrogen, and oxy- 
gen abundances in the solar neighbourhood constrains the 
shape of the IMF (see also Gibson & Mould 1997). Cou- 
pling chemodynamical models to spectrophotometric codes 
(see § 6) allows one to constrain, for example, the zero point 
of the Tully-Fisher relation (e.g. Steinmetz & Navarro 1999; 
Navarro & Steinmetz 2000; Koda, Sofue & Wada 2000; 
Abadi et al. 2003a). 

We now review the methodology of chemodynamical 
codes and their application in modeling the evolution of the 
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Milky Way. The discussion is framed primarily within the 
context of our N-body/Smoothed Particle Hydrodynamics 
software package GCD+ (Kawata 1999; Kawata & Gibson 
2003a,b; Brook et al. 2003a,b). 

5.1 Brief Introduction to GCD+ 

GCD+ was originally (Kawata 1999; Kawata & Gibson 2003a) 
based upon Katz et al.'s (1996) TreeSPH code, combining 
the tree algorithm (Barnes & Hut 1986) for the computa- 
tion of the gravitational forces, with the smoothed particle 
hydrodynamics (SPH) (Lucy 1977; Gingold & Monaghan 
1977) approach to numerical hydrodynamics. The dynam- 
ics of the dark matter (DM) and stars is calculated by the 
N-body scheme, and the gas component is modeled using 
SPH. It is fully Lagrangian, three-dimensional, and highly 
adaptive in space and time owing to individual smoothing 
lengths and individual time steps. Moreover, it includes 
self-consistently important physical processes, such as self- 
gravity, hydrodynamics, radiative cooling, star formation, 
SNe feedback, and metal enrichment. 

Metallicity-dependent radiative cooling is taken into ac- 
count, following Sutherland & Dopita (1993). The cooling 
rate for a gas with solar metallicity is larger than that for 
gas of primordial composition by more than an order of 
magnitude. As such, the cooling provided by metals should 
not be ignored in numerical simulations of galaxy formation 
(Kay et al. 2000). 

GCD+ takes into account the energy feedback and metal 
enrichment from both SNe II and SNe la, as well as en- 
richment from ILMS. We assume that each massive star (> 
8 Mq) explodes as a Type II supernova. The SNe la rates 
are calculated using the model proposed by Kobayashi, 
Tsujimoto & Nomoto (2000). The yields of SNe II, SNe la, 
and ILMS are taken from Woosley & Weaver (1995), Iwamoto 
et al. (1999), and van den Hoek & Groenewegen (1997), 
respectively. The mass, energy, and heavy elements are 
smoothed over the neighbouring gas particles using the SPH 
smoothing algorithm. For example, when the i-th star par- 
ticle ejects a mass of Msn.s, a mass increment is applied to 
the j-th neighbour gas particle, as such: 

AM SN , 3 = ^-M SNli W(r«/hi), (3) 

Ps,i 

where 

Pg,i = (Psim)) ="y]rn j W(r ij /hi), (4) 

and W(x) is an SPH kernel (Kawata & Gibson 2003a). GCD+ 
monitors the evolution of all relevant chemical elements self- 
consistently with the dynamics of the gas, stars, and dark 
matter. The photometric evolution of the stellar popula- 
tions is also calculated simultaneously via the Kodama & 
Arimoto (1997) spectral synthesis package. 

5.2 Disk Galaxy Formation in ACDM Cosmolo- 
gies 

We present now a sample disk galaxy formation simula- 
tion, undertaken within a ACDM cosmological framework, 
in order to demonstrate the capabilities of chemodynam- 
ical codes such as GCD+. We carried out a series of high- 
resolution simulations within a A-dominated cold dark mat- 
ter (ACDM) cosmology (fi =0.3, A =0.7, n b =0.019h~ 2 , 



h=0.7, <7g=0.9). We used a multi-resolution technique (de- 
scribed below) to achieve high-resolution in the region of in- 
terest, including the tidal forces from large scale structures. 
The multi-resolution initial conditions were constructed us- 
ing Bertschinger's (2001) publically available software GRAFIC2. 
First, a low-resolution N-body only simulation of a comov- 
ing 20h~ 1 Mpc diameter sphere was performed; the mean 
separation of the particles therein was 20/i _1 /64 Mpc. The 
mass of each particle was 3.63 x 10 9 M©, and a fixed soft- 
ening of 18.0 kpc was applied. Next, at redshift z — 0, we 
selected an ~10 Mpc diameter spherical region which con- 
tains several galaxy-sized DM halos. We traced the parti- 
cles which fall into the selected region back to the initial 
conditions at z — 43.5 and identified the volume which 
consists of those particles. Within this arbitrarily shaped 
volume, we replaced the low resolution particles with par- 
ticles a factor of 64 times less massive. The initial den- 
sity and velocities for the less massive particles are self- 
consistently calculated by GRAFIC2, taking into account 
the density fields of the lower resolution region. Finally, we 
re-simulated the full volume (20h~ 1 Mpc sphere), but now 
including all gas dynamics, cooling, and star formation. 8 
The surrounding low-resolution region contributes to the 
high-resolution region only through gravity. The mass and 
softening length of individual gas (dark matter) particles in 
the high-resolution region were 7.33 x 10 6 (4.94 x 10 7 ) M 
and 1.14 (2.15) kpc, respectively. 

At z = 0, using a friends-of- friends methodology, we 
identified six stellar systems which consisted of more than 
2000 star particles. Two of these systems resembled large 
disk-like systems with kinematics consistent of rotational 
support; one was chosen as the target galaxy. The total 
virial mass of this galaxy is ~2.4x 10 12 Mq, and there exists 
a companion with ~30% the mass of the target at a galacto- 
centric distance of ~270 kpc. The virial mass is defined as 
the mass within the virial radius, which itself is the radius 
of a sphere containing a mean density of 178 Qq' 45 times 
the critical values (p cr it = 3Ho/87rG), after Eke, Navarro & 
Frenk (1998). 

Figure 5 shows the morphological evolution of the dark 
matter in a central portion of the simulation volume, and 
the evolution of the stellar component in a 200 kpc region 
centred on the target galaxy. The lower panels correspond 
to the predicted J-band (in the rest frame) image of the 
target galaxy. In our simulations, the star particles each 
carry their own age and metallicity "tag", due to the self- 
consistent chemodynamical nature of the calculation. This 
enables us to generate an optical-to-near infrared spectral 
energy distribution for the target galaxy (here, using the 
Kodama & Arimoto 1997 spectral synthesis package). 

The spectral energy distribution of each star particle is 
assumed to be that of a simple stellar population - i.e. a 
coeval and chemically homogeneous assembly of stars. We 
take into account the fc-correction, but do not consider here 
the effects of dust absorption. Figure 5 demonstrates that 
the galaxy forms through conventional hierarchical cluster- 
ing before 2 = 1; the disk has subsequently been built-up 
smoothly. 

Figure 6 demonstrates the time evolution of the phase 
space information of stars within a satellite which accretes 
onto the central galaxy at z ~ 0.3. At z = 0.34, the 

8 The gas component though was included only within the 
high-resolution region. 
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Figure 5: Dark matter density map of a portion of the 17 Mpc (comoving) simulation volume (upper panels), and 
predicted J-band (AB magnitude) image (physical scale) of the target galaxy (lower panels), over the redshift 
range 2=1.55 to z—0. The projection in the lower panels has been chosen in order to view the target galaxy 
edge-on at z=0. 



tidal stream of the satellite appears. The stream is identi- 
fiable in the radial velocity (VR, ga i) - galactocentric radius 
(-Rgai) diagram, in addition to the VR lga i - rotational veloc- 
ity (V^gai) diagram. This satellite is completely disrupted 
at z = 0, and the VR, g ai-i?gai information of the member 
star particles of the satellite also disappear. However, in 
the VR.gai-V^gai diagram, the distribution of the member 
star particles is similar to that at z = 0.34. We confirmed 
that the velocity phase space information of the member 
stars is conserved very well, even after the satellite is spa- 
tially disrupted (see also Helmi & White 1999; Helmi & de 
Zeeuw 2000). 

It is worth noting however that it is difficult to identify 
the member stars of the accreted satellites in the velocity 
phase space diagram of the full sample of observed stars. 
The red dots in Figure 6 represent the halo stars. Here, we 
define the star particles with [Fe/H]< —0.6 (Chiba & Beers 
2000) as the halo stars. Figure 6 shows that in the velocity 
phase space diagram, the distribution of the satellite mem- 
ber stars overlaps that of the halo stars. Thus, to identify 
the member stars of the accreted satellites, additional ob- 
servational information is required. Since the member stars 
formed within a small galaxy, it is anticipated that they 
might each inherit a unique chemical "fingerprint" (Free- 
man & Bland-Hawthorn 2002). Hence, the combination of 
such chemical tags and kinematics can be a powerful tool to 
identify the field stars which orginated within now accreted 
satellites. Comparing such observations and the results of 



chemodynamical simulations will be critical to understand- 
ing the formation history of the Milky Way. Unfortunately, 
current numerical simulations still struggle to overcome the 
classical "over-cooling problem" (White & Frenk 1991), the 
signatures of which are a high-redshift star formation rate 
in excess of that observed, and stellar halos which are both 
too massive and too metal-rich (Brook et al. 2003a,b; Helmi 
et al. 2003). The exact physical mechanism required to 
solve this problem remains uncertain, although a framework 
predicated upon an enhanced supernova feedback efficiency 
is one leading candidate (Navarro & Steinmetz 2000; Brook 
et al. 2004). 



6 Spectrophotometry 

Determining the epoch of galaxy formation and understand- 
ing the consequent chemical evolution are amongst the fun- 
damental quests of modern cosmology. An complementary 
approach to addressing the chemical evolution of galaxies 
is via the use of spectrophotometry. One can directly de- 
rive the age and metallicity of a galaxy, and their respective 
gradients therein, by comparing its observational integrated 
colours and/or its spectral line indices with theoretical pre- 
dictions from stellar population synthesis techniques. The 
integrated properties of star clusters around a galaxy can 
also be used to induce their host galaxy's chemical evolu- 
tionary path. 
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Figure 6: Snapshots of the spatial and velocity distributions of the member star particles (blue dots) of a 
satellite which accretes onto the central galaxy at z ~ 0.3. Left panels show the edge-on projection of the 
accreting satellite as well as J-band (AB magnitude) images of the target galaxy (colour contour image, with 
levels as in Figure 5) . Middle and right panels show the radial velocity versus the galactocentric radius (middle) 
and circular velocity (right) of particles with galactocentric radius less than 50 kpc. The red dots show star 
particles with [Fe/H]<— 0.6, representative of the population of halo stars in the central galaxy. 
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Figure 7: The metallicity [Fe/H] of Milky Way globu- 
lar clusters are well-reproduced from this spectral in- 
dex - index plot. The data shown - including Mg2 
and [MgFe] indices - are are taken from Cohen et al. 
(1998) and the metallicities in parentheses are from 
Harris (1996, Feb. 2003 Version). 



Figure 8: The relatively low-reddened Milky Way glob- 
ular clusters [E(B — V) < 0.2] are used to calibrate 
our models in the (B — V) vs. [Fe/H] plane. The 
dashed and solid lines represent ages of 10 Gyr and 12 
Gyr, respectively. Filled (open) circles correspond to 
inner (outer) halo globulars. 



Necessary ingredients for theoretical spectrophotomet- 
ric predictions include stellar evolutionary tracks, isochrones, 
and a corresponding stellar atmosphere library which covers 
a wide range of stellar parameters, such as metallicity, tem- 
perature, and surface gravity. The stellar population syn- 
thesis models presented below, for example, are based upon 
the Y 2 Isochrones 9 (Kim et al. 2002) with [a/Fe]=+0.3, 
coupled to the post-red giant branch stellar evolutionary 
tracks of Yi, Demarque & Kim (1997). The stellar library 
of Lejeune, Cuisinier & Buser (1998) was taken for the con- 
version from theoretical quantities to observable quantities. 

Both age-sensitive and metallicity-sensitive spectropho- 
tometric quantities are initially constructed for a grid of 
simple stellar populations, for a range of age and metallic- 
ity. Composite stellar population spectrophotometric quan- 
tities can then be generated by convolving any given star 
formation history (with the requisite self-consistent treat- 
ment of chemical enrichment) with the grid of simple stellar 
population results (weighted by the number of stars popu- 
lating each stellar evolutionary stage of each simple stellar 
population). The spectrophotometric quantities calculated 
from these composite populations can then be compared 
directly with observational data, and the age and metallic- 
ity of the underlying stellar population extracted. Once we 
have a fair selection of sample galaxies in terms of age and 
metallicity, we are able to investigate the detailed abun- 
dance properties along the age sequence to understand the 
chemical evolution of galaxies. 

Figure 7 demonstrates the use of spectral index vs. in- 

9 http : // csaweb.yonsei . ac .kr/~kim/yyiso .html 



dex plots as tools for estimating metallicity. A sample of 
integrated spectra for 12 Milky Way globular clusters is 
shown (from Cohen et al. 1998), which was used to de- 
fine a grid of Lick indices which are then compared with 
our recent models (H.-c. Lee & B.K. Gibson, in prepa- 
ration). The [MgFe] index is defined as ^/Mgbx < Fe >, 
where < Fe > is (Fe5270 + Fe5335)/2. It is found that 
the metallicity [Fe/H] that is acquired independently from 
the Harris compilation (Feb. 2003 Version) is surprisingly 
well recovered. This kind of calibration is necessary and 
should be quite useful for the derivation of metallicity for 
extragalactic star clusters and/or external galaxies that are 
not resolved into individual stars. 

Figure 8 shows the importance of the realistic man- 
ifestation of horizontal-branch (HB) morphologies in the 
stellar population synthesis models for relatively old stellar 
systems (r > 8 Gyr) (Lee, Yoon & Lee 2000; Lee, Lee & 
Gibson 2002). It appears that models with a proper treat- 
ment of blue HB stars reproduce the differences between 
inner and outer halo clusters, 10 in the sense that the in- 
ner halo clusters are not only more tightly grouped along 
the isochrone than the more scattered outer halo counter- 
parts, but also relatively older. This is interesting in that 
the outer halo clusters with the wider range of colours at a 
given metallicity may indicate a different chemical "origin" , 
perhaps from different environments such as satellite dwarf 
galaxies. 

Another important parameter that controls spectropho- 

10 The inclusion of blue HB stars also led Lee et al. (2000,2002) 
to suggest that giant elliptical galaxies may be ~ 1—3 billion years 
older than the Milky Way. 
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Figure 9: Our simple stellar population predictions for 
(B — R) vs. M/Lr as a function of age for three IMF 
slopes. The solid line with filled circles corresponds to 
the standard Salpeter (1955) IMF (by number) with an 
exponent %=2.35. The dotted lines with open circles 
and open squares are x=3.35 and x=3.85 cases, respec- 
tively. At a given age and IMF exponent, the metallic- 
ity [Fe/H] is -2.51, -1.90, -1.51, -0.90, -0.49, -0.17, 
+0.17, and +0.39, from left to right. The filled pen- 
tagons are Fuchs (2002) sample of LSB galaxies. From 
the left panel, it is suggested that these LSB galaxies 
with high M/Lr ratios and blue colours ((B — R) < 
1.0) may be consistent with a relatively low metallic- 
ity and rather recent (<2 Gyr ago) star formation, by 
introducing a steep IMF (Lee et al. 2004). 



tometric quantities, as well as the mass-to-luminosity (M/L) 
ratio, is the IMF. A top-heavy IMF leads to highly efficient 
chemical enrichment due to a preponderance of Type II 
SNe. Conversely, chemical enrichment is minimised under 
the adoption of a bottom-heavy IMF. A recent attempt by 
Fuchs (2002) to set constraints on the mass of the disks 
of low surface brightness (LSB) galaxies by employing den- 
sity wave theory is intriguing in this respect. Fuchs used a 
sample of five LSBs with clear spiral structure to claim that 
each possess surprisingly high stellar mass-to-luminosity ra- 
tios in the R-band (M/Lr > 3), in addition to their blue 
colours ((B — R) < 1.0). Figure 9 demonstrates that the 
range of colours and mass-to-luminosity ratios spanned by 
the Fuchs LSB sample are consistent with a low-metallicity 
([Fc/H]< — 1.5), recent (<2Gyr), burst of star formation, 
under the assumption of an IMF significantly steeper than 
that of Salpeter (1955). 



7 Future Directions 

Any model for the formation and evolution of the Milky 
Way is only as good as the observational data upon which 
it is calibrated. We wish to end this review with an outline 
of some of the exciting observational programs which will 
come to fruition over the coming decade. Each of these 
datasets is capable of constraining - in a new and significant 
manner - the GCE models discussed in § 3 — § 6. 

We are fortunate that for the Milky Way we can, in 
principle, obtain full six-dimensional phase space (spatial 
and kinematical) and chemical information, for individual 
stars. Such data provides unique insights into the detailed 
formation history of our Galaxy (e.g. Eggen, Lynden-Bell & 
Sandage 1962; Chiba & Beers 2000; Bekki & Chiba 2001; 
Freeman & Bland-Hawthorn 2002; Brook et al. 2003a). 
The discovery of the disrupting Sagittarius dwarf (Ibata 
et al. 1994) and other halo substructure signatures (Helmi 
et al. 1999; Chiba & Beers 2000; Gilmore, Wyse & Nor- 
ris 2002; Brook et al. 2003a) has also demonstrated the 
value of stellar kinematics in reconstructing satellite accre- 
tion events. Such events can (potentially) aid in determin- 
ing the fraction of the halo which was accreted and the 
fraction which formed in situ (Helmi & White 1999; Helmi 
& de Zeeuw 2000; Harding et al. 2001). 

We identify three landmark projects which will advance 
significantly the field of "Galactic Archaeology" : 

• RAVE 11 (RAdial Velocity Experiment) is an 

ambitious all-sky survey (complete to V=16) aimed 
at measuring the radial velocities (with precision <2 km s _1 ) 
metallicities, and abundance ratios (both to ~0.1 dex 
precision) of 50 million stars using the United King- 
dom Schmidt Telescope (UKST), together with a north- 
ern counterpart, over the period 2006 - 2010, provid- 
ing a vast stellar kinematic and chemical database. 
A two-year pilot survey commenced on the UKSL 
in April 2003, making use of the existing 6dF Facil- 
ity. This pilot survey will obtain comparable qual- 
ity data to the main survey for 10 stars (spanning 
~8500deg 2 , 9<I<12), of which approximately half 
have accurate Tycho-2 proper motions. 



http: //astronomy . swin. edu . au/RAVE/ 
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• ESA's GAIA satellite mission, scheduled for 
launch in 2010, will conduct a census of ~10 9 stars - 
~100 revisits per star over a five-year period - mea- 
suring positions of all objects down to V=20. Posi- 
tional accuracies of ~4/xas at V=10 and ~10pas at 
V=15 will yield distances accurate to 10% at helio- 
centric distances of ~10kpc. Radial velocity accu- 
racy comparable to that of RAVE is expected (also 
to V=16). 

• NASA's Space Interferometry Mission - SIM 13 . 

scheduled for launch in 2009, complements GAIA's 
4 /xas positional accuracy with a narrow-angle mode 
allowing accuracies of ~ 1 pias. 

There currently exists little compelling evidence that 
a massive (>10 10 Mq) satellite has been accreted by the 
Milky Way within the past ~10 Gyr (Gilmore et al. 2002) . 14 
Instead, later accreted satelites are more likely to resemble 
dwarf systems such as the Local Group's dwarf spheroidals 
(dSphs). Recent (and spectacular) high- resolution spec- 
troscopic observations of individual stars in dSphs have 
demonstrated that their ratio of a-elements to iron tend 
to be close to solar (e.g. Shetrone et al. 2001,2003). This 
elemental abundance pattern is not reflected in the present- 
day halo field stars - a somewhat surprising result consid- 
ering the aforementioned kinematical evidence for halo ac- 
ccretion events. An inescapable conclusion to be drawn 
from these data is that the present-day Local Group dSphs 
are not the primary "stellar donors" to the Galactic halo 
(Tolstoy et al. 2003). Whether such accreted satellites are 
responsible for the thick disk though remains an intrigu- 
ing possibility (Bekki & Chiba 2002; Abadi et al. 2003b). 
Chemical "tagging" of Local Group dSphs is still very much 
in its infancy, but a burgeoning field which will see enor- 
mous increases in sample sizes over the coming few years, 
thanks to substantial investments of 8m-class time. 

The hunt for the most primitive stars in the Milky Way 
has led to some of the most exciting discoveries in recent 
astronomy, including the detection of a star with an iron 
abundance less than 1/200,000 that of the Sun (Christlieb 
et al. 2002). Having formed from almost pristine primordial 
gas, this star (HE0107-5240) is a nearby counterpart to the 
high-redshift universe and provides insight into the earliest 
epochs of Galaxy formation. Two factors could prohibit 
the detection of a bona fide zero- met allicity Pop III star: 
1) metal-free gas might favour the formation of higher mass 
stars with short lifetimes, such that there are no surviving 
Pop III stars today (Nakamura & Umemura 2001); and 2) 
the surfaces of Pop III stars may have been polluted with 
metals either from internal processing or through the accre- 
tion of interstellar matter (Shigeyama et al. 2003). Despite 
their paucity, thousands of candidate metal-poor stars have 
been selected for follow-up spectroscopy by programs such 
as the HK and Hamburg/ESO surveys (Beers 2000). The 
tally of ~100 stars with [Fe/H]<-3 found by the HK sur- 
vey is expected to grow by more than a factor of five with 
the Hamburg/ESO program (Christlieb 2003). As alluded 
to in § 4, the abundance pattern in extremely metal-poor 

12 http : //astro . estec . esa.nl/GAIA/ 
13 http : //sim. j pi .nasa.gov/ 

14 The potential impact that such a massive satellite would 
have upon the disk supports this contention (Kawata et al. 
2003). 



stars may reflect the chemical fingerprint of a single Pop III 
star and provides empirical constraints on models of "The 
First Stars" (e.g. Umeda & Nomoto 2003). A project us- 
ing the Sloan Digital Sky Survey dataset to find metal-poor 
stars in the halo and thick disk of the Milky Way is also 
underway (Allende Prieto et al. 2003). 

We have only touched the "tip of the iceberg" here - 
the combination of "big glass" and "big surveys" over the 
next decade will see an explosion of chemo-kinematical data 
coming on-line. lj Such a wealth of impending information 
makes this the optimal time to be a theorist interested in 
understanding the detailed formation and evolution of our 
Milky Way Galaxy. 
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